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Abstract: This study introduces a new signal analysis method called SCSA, based 
^ on a semi-classical approach. The main idea in the SCSA is to interpret a pulse- 

shaped signal as a potential of a Schrodinger operator and then to use the discrete 
spectrum of this operator for the analysis of the signal. We present some numerical 
examples and the first results obtained with this method on the analysis of arterial 
blood pressure waveforms. 
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r^ 1 Introduction 

o 
o . 

This paper considers a new signal analysis method where the main idea is to 
interpret a signal y as a multiplication operator, <p — ¥ y.(p, on some function space. 
The spectrum of a (formally) regularized version of this operator, denoted Hh{y) and 
defined by 



X 



H h (y)ifj = -h 2 —^- yi p, iPEH 2 (R), h>0, (1) 

ax 
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for small h, is then used for the analysis instead of the Fourier transform of y. 
Here H 2 (R) denotes the Sobolev space of order 2. In this method, the signal is 
interpreted as a potential of a Schrodinger operator. This point of view seems useful 
when associated inverse spectral problem is well posed as it will be the case for some 
pulse shaped signals. 

We define Hh on a space B such that 

B = {yeL\(R), y(x)>0, WxeR, ^eL\R), m = l,2}, (2) 

/+oo 
|V(x)|(l + \x\)dx < oo}. L\(R) is known as the Faddeev 
-oo 

class [9]. 

For A < 0, we denote Nh(X; y) the number of eigenvalues of the operator Hh{y) 
below A. Under hypothesis (J2J), there is a non-zero finite number Nh = Nh(0; y), as 
it is described in proposition 12. II We denote — K 2 nh the negative eigenvalues of Hh(y) 
with n nh > and K lh > K 2h > ■ ■ ■ > K nh , n = 1, • ■ • ,N h . Let ip nh , n = 1, • ■ • , N h be 
the associated L 2 -normalized eigenfunctions. 

In this study, we focus our interest in representing the signal y with the discrete 
spectrum of Hh(y) using the following formula (J3J) 

N h 

y h (x) =4h^K nh if)l h (x), xeR. (3) 

n=l 

As we will see, the parameter h plays an important role in our approach. Indeed, 
as h decreases as the approximation of y by y^ improves. Our method is based on 
semi-classical concepts i.e h — > and we will call it Semi-Classical Signal Analysis 
(SCSA). 

In the next section, we will present some properties of the SCSA. In section 3, we 
will consider a particular case of an exact representation for a fixed h and show its 
relation to a signal representation using the so called reflectionless potentials of the 
Schrodinger operator. Section 4 will deal with some numerical examples and section 
5 will present some results obtained on the analysis of Arterial Blood Pressure (ABP) 
signals using the SCSA. A discussion will summarize the main results and compare 
the SCSA to related studies. In appendices A, B and C some known results on direct 
and inverse scattering transforms are presented. 



2 SCSA properties 

To begin, we focuss our attention on the behavior of the number N^ of negative 
eigenvalues of Hh(y) according to h as described by the following proposition. 

Proposition 2.1. Let y be a real valued function satisfying hypothesis (0|). Then, 

i) The number Nh of negative eigenvalues of Hh(y) is a decreasing function of h. 

ii) Moreover if y E L^iM.), then, 

I /-+00 

lim hNh = — \Jy(x)dx, (4) 

h ^° 7T j_ oe> 

Proof. 

i) The proof of this item is based on the following lemma. 

Lemme 2.1. t2j]28j For A < 0, let iVi(A; V) be the number of negative eigen- 
values of Hi(y) less than A. Let V and W be two potentials of the Schrodinger 
operator such that W < V, then 

iVi(A;W)<iVi(A;V), VA < 0. (5) 

Let y e B. We put V = —ry, W = —ry, with < h 1 < h 2 . 

h{ h 2 2 

For A = 0, we obtain iVi(0; -r^y) < ^(0; -r^y), and we have iVi(0; -hy) = N h ., 

/i 2 h^ j 

j = 1,2, which proves the result. 

ii) Let y G B D L^(M.) and A e] — y ma x, 0]. We denote 

S,(h,X)= J2 (A + ^,) 7 , 7>0 (6) 

K 2 , <A 

nh — 

the Riesz means of the values — K^ h ^ ess than A. Remark that So(h,X) = 
N h (X;y). 

Property ii) results from the following lemma 12.21 



Lemme 2.2. JMBj For y G L 7+ 5(R), y(x) > 0, Vx G R and 7 > 0, we ftave 

/+00 
y( X y + Ux, (7) 

-00 

where L C J: is the classical constant given by 

Li 



r( 7 + i; 



"i 2 v ^r( 7 + |)' 

for all 7 > 

By taking 7 = in ((Tj) we get the result. 
Let us now study some properties of the negative eigenvalues —K^ h , n = 1, • • • , iVj, 

ofiZhfo). 

Proposition 2.2. Lei y G C°°(R), zm'tfi y(x) > 0, Vx G R and sitc/i that 3~f G R, 
minR (— y + 70) > and\/a G N, 3C a > such that |^f | < C Q (— y + 7o), i^en every 
regular value of y is an accumulation point of the set (i < ? nh , n = 1, • • • , Nh) (v is a 
regular value if0<v< y max and ify(x) = v then \-^p-\ > 0). 

Proof. 

We want to show that every regular value of y is an accumulation point for the 
set (n^h, n = 1, • • • , iV/J. For this purpose, we use the following result shown by B. 
Helffer and D. Robert [33]. 

Theorem 2.1. \L$ Let y G C°°(K), with y{x) > 0, \/x G R such that for 70 G R, 
miniR (—y + 70) > and for all a E N, there is a constant C a > such that \g^\ < 
C a {—y + 7o)- -^ A < h mm f|x|-^oo( — y{ x )), then for < 7 < 1 the Riesz means (0|) 
are given by 
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+00 



S,(h, A) = ^ ( L* J I A + y(x)|f * dx + 0(/i 1+ ^) J , (9) 

where \V\ + is the positive part ofV and Uz, known as the classical constant, is given 

bym- 



We put 7 = in (J6j). We notice that So(h, A) = Nh(X,y). Substituting 7 by in 
(ED, we get 



1 r + °° 1 

S (h, A) = -L* J y/\X + y(x)\ + dx + 0(1). (10) 

We suppose that there is a regular value y of y that is not an accumulation point 
of the set (K^ h , n = 1, • • • , Nh). So there is a neighborhood V(yo) of y$ and a value 
ho, small enough such that V/i < /io, V(y ) does not contain any element element K^ h . 

Moreover, we can choose V(yo) small enough such that 

inf{|Jjp2| } for x such that n(i) 6 %)} = c > 0. (11) 

Then, we can take V(yo) = [yi,V2[i with y\ and 2/2 some regular values of y and 
< Vi < 2/2- 

For all h < ho, the difference So(h, —yi) — So(h, —y 2 ) represents the number of 
elements of (—K^ h , n — 1, ■ ■ ■ , Nh) in the interval ] — y 2 , —yi]- However this set is 
empty because there is no element in the neighborhood of y . 

Denoting X(X) = {x\y(x) + A > 0}, we have from ( ITU]) 
So(h,-yi)-S (h,-y 2 ) = t l o( Vv( x ) ~ Vidx- / y/y(x) - y 2 dx 

n \Jx(- yi ) Jx(-ya) 

+ 0(1), (12) 

so, as the left quantity is null, we obtain 

VM^i* - I VM^i*, (13) 

X(- W ) JX(-y 2 ) 

and we have X(-y ± ) = X(-y 2 ) U y~ 1 ([yi,y2\), so 



(Vy( x ) -y\- Vy( x ) - vi)dx + / Vy( x ) - Vidx = o. (15) 

'x(-y 2 ) Jy~ 1 ([yi,y2[) 

Hence as these two integrals are positives, we get 

/ ym^* = 0, (16) 

Jy HbiiwD 



therefore \y(x) — yi\ = almost everywhere in y~ l ([yi, 2/2D, then y\ = y 2 , which is a 
contradiction. 

Now, proposition 12.21 introduces an interesting property of the SCSA. Indeed, 
remembering that 



, d 2 ib nh (x) 
dx 2 



-p-2 y(x)ifj nh (x) = -K 2 nh ip nh {x). (17) 



and by multiplying the previous equation by ip n h(x) and integrating it by part we 
get 

K nh = / y^nh( X ) dx - h / (— ^ ^- ( 18 ) 

Then, we notice that —y max < — K^ < as it is illustrated in figure [TJ Hence, for 
a fixed value of h, i^ nh can be interpreted as particular values of y which define a 
new quantization approach that can be interpreted by semi-classical concepts. The 
SCSA appears then as a new way to quantify a signal. 




d 2 
Figure 1: The negative spectrum of the Schrodinger operator —h 2 -— — y(x) provides 

ax 2 
a quantization of the signal y 



To finish this section, we examine the convergence of the first two momentums of 
n n h when h — > through proposition 12.31 These quantities could be very interesting 
in signal analysis as it is mentioned in section 5 and also in a recent study |20j . 



Proposition 2.3. Under hypothesis $B), we have 

N h i /. + OO 



and remembering that 



J ' "_ i r+oo 

\imh^K nh = - y(x)dx, (19) 

"^ 71=1 ^~°° 



y h (x) = 4hJ2^nh^ 2 nh (x), xeR, (20) 

71=1 



we have 

»+oo 



lim / yh(x)dx = / y(x)dx. (21) 

Moreover if y G £ 2 (M) 7 i/ien, 

h ft E^ = ifi/ ^^ ( 22 ) 

^o ^ i6 y_^ 
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Proof. The limits (fT9|) and (j2"2"|) are deduced from lemma [2721 for 7 = - and 7 = - 

respectively, 

\imhy^K nh = - y(x)dx, (23) 

71=1 J ~°° 

^ q /-+00 

iTo h ^2 K3nh= iQj y 2 ( x ) dx - ( 24 ) 



71=1 

By integrating ([3]) we have 



/+00 N h 

y h (x)dx = 4h'Y^K nh . (25) 

-°° n=l 



Then, combining U25J and (H2J, we get (12TJ) . 



3 Exact representation and reflectionless potentials 

In this section, we are interested in an exact representation of a signal for a fixed 
h and its relation to reflectionless potentials (reflectionless potentials are defined 
in the appendix B) of the Schrodinger operator as it is described in the following 
proposition: 

Proposition 3.1. The following properties are equivalent, 

i) Equality in (C0|) holds for a finite h ; 

ii) 3h such that yh = y ; 

V 
Hi) 3h such that — is a reflectionless potential of HuV). 

h 2 

Proof 

The following proof uses some concepts and results from scattering transform 
theory that are recalled in the appendix. First, we suppose that i) is fulfilled then 

/ + OO N h 

y{x)dx = Ah^K nh . (26) 

■°° n=\ 

y 

Writing the first invariant [56] (see the appendix C) for the potential -— we have 

AT 



N h 



(27) 



/+oo "ft i2 r+oo 
y(x)dx = 4h^K nh -\ / In (1 - \R r{ i )h (k)\ 2 )dk, 
■oo n=1 7T J-oo 

where R r ^h{k) is the reflection coefficient (see the appendix A). 
Using (125]) . we get 

/+oo r+oo r 2 p+oo 

y{x)dx = / y h {x)dx + — / In (1 - \R r(l)h {k)\ 2 )dk. (28) 

-oo J— oo " J~ oo 

Then, from ([30]) and (EH]) we obtain 

/+oo 
ln(l-\R r{l)h (k)\ 2 )dk = 0. (29) 

-oo 

The reflection coefficient of a Schrodinger operator satisfies \Rr(i)h{k)\ < 1, Vfc G I 
(see for exemple [6]). Then, we get ln(l — \R r ^h{k)\ 2 ) < 0. Equality ([29]) is then 



fulfilled if and only if ln(l — \R r ^h{k)\ 2 ) = 0, k G M. a.e; which is true only if 
\Rr(i)h{k)\ = 0, k G M a.e. This property defines a reflectionless potential. So, 
z) =>- Hi). 

Now, using the Deift-Trubowitz formula (152|) (see appendix B) that we rewrite 

for the potential — — and taking R r ^h(k) = 0, we can deduce that statement (iii) 

implies statement (ii). 

Then, if we suppose that i/h = y for a given value of h we have 



Vh(x)dx = / y(x)dx, (30) 

J — oo J —oo 

hence ii) =>- z) 

4 Numerical results 

In this section, we are interested in the validation of the SCSA through some 
numerical examples. For this purpose, it will be more convenient to consider the 

y 

problem associated to HA—). Therefore, in order to simplify the notations, we put 

h 2 

1 k 2 

— = X, N h = N x and -^ = k£ y , n = 1, • • ■ ,N X . We denote the L 2 -normalized 
hr h 2 

eigenf unctions ^ nx , n — 1, • • • , N x . Formula ([3]) is then rewritten 

4 Nx 

y ^ = :E «w^n X ( x )> x G R ' ( 31 ) 

^ n=l 

We start by giving the numerical scheme used to estimate a signal with the 
SCSA. Then, the sech-squared function will be considered. This example illustrates 
the influence of the parameter \ on the approximation. Gaussian, sinusoidal and 
chirp signals will be also considered. 

4.1 The numerical scheme 

The first step in the SCSA is to solve the spectral problem of a one dimensional 
Schrodinger operator. Its discretization leads to an eigenvalue problem of a matrix. 
In this work, we propose to use a Fourier pseudo-spectral method pj)], |30]. The 
latter is well-adapted for periodic problems but in practice it gives good results for 
some non-periodic problems for instance, rapid decreasing signals. 



We consider a grid of M equidistant points Xj, j — 1, • ■ ■ , M such that 

a = x± < %2 < • ■ ■ < xm-i < %m = b. (32) 

Let Ax = be the distance between two consecutive points. We denote y~ and 

M - 1 3 

ipj the values of y and ip at the grid points Xj, j — 1, ■ • • , M 

Vj = V( x j)i i>i = ^(^i). 3 = li • • • » M - ( 33 ) 

Therefore, the discretization of the Schrodinger eigenvalue problem leads to the 
following eigenvalue matrix problem 

{-D 2 - X &ag{Y))± = \±, (34) 

where diag(V) is a diagonal matrix whose elements are yj, j = 1, • • ■ ,M and tp = 
["01 '02 5 ••• "0M-1 V'm] • -D2 is the second order differentiation matrix given by [30 j, 
• If M is even 

f=^-i for*-; 

A 2 3A 2 6 ~ J ' 

Da(fc ' J) = 7A^ - (-l) fc - 1 * . for k 4 7. (35) 



2 • 2 ( (k-j)A 

sin ' A — ="— 



If M is odd 



7T 2 1 



3A 2 12 P ° Ur k = h 

~~W?\ -(-l) fc "^ * .^ cotfe^) pour Mi, 



sm 



2 D ;„ r ( fc -j) A 



(36) 

2tt 
with A = — . the matrix D 2 is symmetric and definite negative. To solve the 

eigenvalue problem of the matrix (— D 2 — xdiag(F)) we use the Matlab routine eig. 
The final step in the SCSA algorithm is to find an optimal value of the parameter 
X- So we look for a value \ that gives a good approximation of y with a small number 
of negative eigenvalues. From the numerical tests, we noticed that the number N x 
is in general a step by step function of \- So, we optimize the following criteria in 
each interval [xi,X2] where N x is constant, 

lA, ,, 4 



J M = m 53 G/< ~ y x^ 2 > Vxi = ~^2 K n X ^li X , i=l,---,M. (37) 

10 



i=l * n=l 
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Figure 2: The SCSA algorithm 



For x large enough (equivalently h small enough), we know an approximate relation 
between the number of negative eigenvalues and x thanks to proposition 12.11 Then, 
we can deduce approximate values of Xi an d X2 according to a given number of 
negative eigenvalues. Fig. [2] summarizes the SCSA algorithm. 

Remark: In practice, we often omit the optimization step and just fix N x to a 
large enough value. 



4.2 The sech-squared function 

In order to illustrate the influence of the parameter x on the SCSA, we first study 
a sech-squared function given by 



y(x) = sech 2 (x — x ), x £ M. 



(38) 



11 



The potential of the Schrodinger operator Hi(xy) is given in this case by: —\ sech [x— 
Xq). This potential is called in quantum physics Poschl-Teller potential. 

It is well-known that the Poschl-Teller potential belongs to the class of reflec- 
tionless potentials if, 

x = Xp = N(N + l), # = 1,2,3,..., (39) 

N being the number of negative eigenvalues of Hi(xy) |22] . 

So, for example, if x — 2, the Schrodinger operator spectrum is negative and 
consists of a single negative eigenvalue given by A = — 1. If x = 6, there are two 
negative eigenvalues: Ai = —4, A2 = — 1 and so on. 

Let us now apply the SCSA to reconstruct y. For this purpose, we must truncate 
the signal and consider it on a finite interval so that the numerical computations 
could be possible. 

Figure [3ja illustrates the variation of the mean square error and N x according to 
X- We notice that, 

• N x is an increasing function of x as described in proposition (12. 1ft . Moreover, 
N x is a step by step function. 



There are some particular values of x f° r which the error is minimal. These 
values are in fact the particular values Xp = N X (N X + 1), N x = 1,2, •• ■ for 
which — XV is a reflectionless potential. 

For all e > 0, there is a value x = Xt such that Vx > Xe> J(x) < e - 



Figure [3jb illustrates the variation of the first four eigenvalues of the matrix 
— D2 — xdiag(F), settled in an increasing way, according to x- We notice that these 
eigenvalues, initially positive, are decreasing functions of x and at every passage from 
N x to N x + 1, a positive eigenvalue becomes negative. 

Otherwise, in figure HJ the first four squared eigenf unctions ip^ , n = 1, • • • , 4 are 
represented for x = 20. Each ^ has n — 1 zeros. 

Figure [5] shows a satisfactory reconstruction of y for N x = 1, 2, 3 and 4. 
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(a) (b) 

Figure 3: (a) Mean square error and number of negative eigenvalues according to 
X for y(x) = sech 2 (x — 6) in [0,15]. (b) Four first eigenvalues according to x f° r 
y(x) = sech 2 (x — 6) in [0, 15] 





Figure 4: Four first eigenfunctions if) 2 for x = 20 for y(x) = sech (x — 6) in [0, 15] 
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Figure 5: Estimation of y(x) = sech (x — 6) in [0,15], from the left to the right: 



N X = 1,N X = 2 (top). N x = 3, N x = 4 (bottom) 
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4.3 Estimation of some signals 



In this section, we are interested in the estimation of some signals with the SCSA. 
In each case, we represent the estimation error, the number of negative eigenvalues 
according to \ an d the real and estimated signals for different values of \. 

We start with a gaussian signal given by: 



1 



(x-nV 



y(x) 



a' 



^k 



For the numerical tests we take a — 0.1 and ^ = 0.75. 

Figures [6] and [7] illustrate the results. We notice that with N x 
is satisfactory and as N x increases better is the approximation. 



(40) 



2, the estimation 
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Figure 6: Mean square error and number of negative eigenvalues according to x f°r 
a gaussian signal 



Now we are interested in a sinusoidal signal defined in a finite interval /, 

( \_ [ A$m(ux + 4>) xel , s 



This signal has negative values, so to apply the SCSA, we must translate the signal 
by Vmin = ~A such that y — y m i n > 0. The Schrodinger operator potential to be 
considered is then given by — x{y — ymin)- F° r the numerical tests, we took A = 2, 
u = n and 6 = —0.5. 



15 






Figure 7: Estimation of a gaussian signal, from the left to the right: N x = 1, N x = 2 
(top). iV x = 3, N x = 4 (bottom) 
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The results are represented in figures [HJ [9] and [101 In [9J a single period of the 
signal is considered while in [101 four periods are represented. In the last case we 
noticed that the negative eigenvalues are of multiplicity 4, they are repeated in each 
period. 



0.02- 
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Figure 8: Mean square error and number of negative eigenvalues according to x f° r 
a sinusoidal signal 



Finally, figures [TT] and [12] illustrate the results obtained in the case of a chirp 
signal. We recall that a chirp signal is usually defined by a time varying frequency 
sinusoid. In our tests, we considered a linear variation of the frequency. 

5 Arterial blood pressure analysis with the SCSA 

ABP plays an important role in the cardiovascular system. So many studies were 
done aiming at proposing mathematical models in order to understand the cardio- 
vascular system both in healthy and pathological cases. Despite the large number 
of ABP models, the interpretation of ABP in clinical practice is often restricted to 
the interpretation of the maximal and the minimal values called respectively the 
systolic pressure and the diastolic pressure. None information on the instantaneous 
variability of the pressure is given in this case. However, pertinent information can 
be extracted from ABP waveform. The SCSA seems to provide a new tool for the 
analysis of ABP waveform. This section presents some obtained results. 
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Figure 9: Estimation of a sinusoidal signal, from the left to the right: N x = 4, N x = 6 
(top). N x = 8, N x = 10 (bottom) 




Figure 10: Estimation of 4 sinusoidal signal periods with N x = 10 
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Figure 11: Mean square error and number of negative eigenvalues according to \ for 
a chirp signal 







Figure 12: Estimation of a chirp signal, from the left to the right: N x = 30, N x = 50 
(top). N x = 80, N x = 100 (bottom) 
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We denote the ABP signal P and P its estimation using the SCSA such that 

N x 

£(*) = -£«■*<*(*)> (42) 

* n=l 

where — k^ , n — 1, • • • , N x are the N x negative eigenvalues of Hi(xP) and ip nx the 
associated L 2 — normalized eigenf unctions. As ABP signal is a function of time, we 
use the time variable t in the Schrodinger equation instead of the space variable x. 
The ABP signal was estimated for several values of the parameter \ and hence 
N x . Figure [T3l illustrates the measured and estimated pressures for one beat of ABP 
and the estimated error with N x = 9. Signals measured at the aorta level and the 
finger respectively were considered. We point out that 5 to 9 negative eigenvalues 
are sufficient for a good estimation of ABP signals [T7], [19J. 

A first interest in using the SCSA for ABP analysis is to decompose the signal 
into its systolic and diastolic parts respectively. This application was inspired from 
a reduced model of ABP based on solitons introduced in |3], [18]. Solitons are in fact 
solutions of some nonlinear partial derivative equations for instance the Korteweg-de 
Vries (KdV) equation which was considered in this reduced model [I]. This model 
proposes to write the ABP as the sum of two terms: and N-soliton, solution of the 
KdV equation describing fast phenomena that predominate during the systolic phase 
and a two-elements windkessel model that describes slow phenomena during the 
diastolic phase. Moreover the KdV equation can be solved with the Inverse Scattering 
Transform (1ST) whose definition is recalled in appendix A. In this approach, the 
KdV equation is associated to a one dimensional Schrodinger potential parameterized 
by time where the potential is given by the solution of the KdV equation at a given 
time. Therefore a relation between the Schrodinger operator and solitons was found 



|10| : solitons are reflectionless potentials. Then according to proposition 13. 1[ the 
SCSA coincides with a soliton representation of a signal for a finite \ when — XV is 
an A^-soliton, where N x denotes the number of negative eigenvalues l of Hi(xy)- 
So each spectral component represents a single soliton. We know that solitons are 
characterized by their velocity which is determined by the negative eigenvalues — KjL, 
n — 1, ■ ■ ■ , iV x of the Schrodinger operator. The largest values K nx characterize fast 
components and the small values of k„ x characterize slow components. From these 
remarks, we propose to decompose equation ( 13"T|) into two partial sums: the first one, 
composed of the N s (N s = 1, 2, 3 in general) largest k„ x and the second partial sum 
composed of the remaining components. The first partial sum describes the systolic 



^ach soliton is characterized by a negative eigenvalue 
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2 22 2.4 2.6 




fa) Aorta 




2 22 2.4 2.6 




(b) Finger 

Figure 13: Reconstruction of the pressure at the aorta and the finger level with the 
SCSA and N x = 9. On the left the estimated and measured pressures, On the right 
the relative error 
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phase and the second one describes the diastolic phase. We denote P s and Pd the 
systolic pressure and the diastolic pressure respectively estimated with the SCSA. 
Then we have 



N e 



N, 



£*(*) = -I>nxVi(t), 
X n=l 



Pd ^ = ~ J2 K n X 1p 2 nx (t). 



X 



(43) 



n=N s +l 



Figure [TH represents the measured pressure and the estimated systolic and di- 
astolic pressures respectively. We notice that P s and Pd are respectively localized 
during the systole and the diastole. 




Estimated diastolic pressure 
Measured pressure 




2.2 2.4 



Figure 14: On the left the estimated systolic pressure. On the right estimated 
diastolic pressure 



6 Discussion 



The spectral analysis of the Schrodinger operator introduces two inverse prob- 
lems: an inverse spectral problem and an inverse scattering problem. 

On the one hand, the inverse spectral problem aims at reconstructing the po- 
tential of a Schrodinger operator with its spectrum (spectral function). It has been 
extensively studied for instance by Borg, Gel'Fand, Levitan and Marchenko [TT| or 
more recently by Ramm [27] . They considered the half-line case and used two spectra 
of the Schrodinger operator with two different boundary conditions in order to recon- 
struct the potential. The inverse spectral problem for a semi-classical Schrodinger 
operator Hh{y) have been recently considered for example by Colin de Verdiere [5] 
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who proposed to reconstruct the potential locally with a single spectrum or Guillemin 
and Uribe p3] who showed that under some assumptions, the low-lying eigenvalues 
of the operator determine the Taylor series of the potential at the minimum. 

In this work we have studied an inverse spectral approach that is different from 
classical inverse spectral problems. Indeed, we used more information to reconstruct 
the potential by including the eigenfunctions as illustrated by equation ([3]). 

On the other hand, the inverse scattering problem aims at recovering the potential 
from the scattering data (see appendix A). Many studies considered this question 
for instance those of Marchenko [25j who proved that under some conditions on 
the scattering data, a potential in ^(R) can be reconstructed from these scattering 
data and gave an algorithm for recovering the potential. We can also quote works 
of Faddeev [9], Deift and Trubowitz [6] considering potentials in Lj(IR), Dubrovin, 
Matveev and Novikov jT] for periodic potentials. 

The convergence of the SCSA when h —)■ is not easy to study. Using the Deift- 
Trubowitz formula (152|) we have tried to consider this problem. However, despite 
interesting results obtained regarding the convergence of some quantities depending 
only on the continuous spectrum of the Schrodinger operator to zero, we did not 
succeed in finding a result of convergence for y^. A study supported by semi-classical 
concepts is now under consideration. The latter is based on the generalization of the 
results of G. Karadzhov |16) . 

The SCSA method has given promising results in the analysis of arterial blood 
pressure waveforms. More than a good reconstruction of these signals, the SCSA 
introduces interesting parameters that give relevant physiological information. These 
parameters are the negative eigenvalues and the invariants which are given by the 
momentums of K nx , n = !,■■■ ,N X introduced in proposition 12.31 For example, 
these new cardiovascular indices allow the discrimination between healthy patients 
and heart failure subjects [21] . A recent study [20] shows that the first systolic 
momentum (associated to the systolic phase) gives information on stroke volume 
variation, a physiological parameter of great interest. 

We have seen in proposition 13.11 that for a fixed value of h the SCSA coincides 
with a reflectionless potentials approximation. This point seems to be an interesting 
avenue of research. Indeed, thanks to the relation between reflectionless potentials 
and solitons, the approximation by reflectionless potentials could have interesting 
applications in signal analysis and in particular in data compression. As we said in 
the previous section, solitons are reflectionless potentials of the Schrodinger opera- 
tor. Gardner et al [10] showed that an A r -soliton is completely determined by the 
discrete scattering data and in particular by 2N parameters which are the negative 
eigenvalues and the normalizing constants. Hence, if — xy is an A^-soliton, it is given 
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by the following formula: 



y{x) 



2JP 

X dx- 



ln(det(J + A x ))(x), xe 



A x (x) is an N x x N x matrix of coefficients 



A x (x) 



Cmx C nx 



\ ^ ?n x~i~^ n X t ^ 



_ n "mx 



+ K- 



"X 



1-1,771 = 1, 



>Nx> 



(44) 



(45) 



where — k 1 and c 



-nx, 



V 



, ■ ■ • , N x are the negative eigenvalues and the normaliz- 
ing constants of Hi{xy)- Hence, this formula provides a parsimonious representation 
of a signal. 

The convergence of the approximation by solitons (or reflectionless potentials) 
when x ~~ >* +°° (equivalently h — > 0) was studied by Lax and Levermore [24] in a 
different context. Indeed they studied the small dispersion limit of the KdV equation 
and approached the initial condition of the KdV equation by an iV-soliton that 
depends on the small dispersion parameter which is in our case h. They showed 
the results in the mono-well potential case and affirmed without prove that the 
result remain still true for multi-well potentials. However the main limitation of this 
approach is the difficulty to compute the normalizing constants. This difficulty can 
be explained by the fact that these constants are defined at infinity which can not 
be handled in the numerical implementation. At the best of our knowledge there is 
no study enabling the computation of these parameters apart a recent attempt by 
Sorine et al 129 1. 



7 Conclusion 

A new method for signal analysis based on a semi-classical approach has been 
proposed in this study: the signal is considered as a potential of a Schrodinger 
operator and then represented using the discrete spectrum of this operator. Some 
spectral parameters are then computed leading to a new approach for signal analysis. 
This study is a first step in the validation of the SCSA. Indeed, we have assessed here 
the ability of the SCSA to reconstruct some signals. We have studied particularly a 
challenging application which is the analysis of the arterial blood pressure waveforms. 
The SCSA introduces a novel approach for arterial blood pressure waveform analysis 
and enables the estimation of relevant physiological parameters. A theoretical study 
is now under consideration regarding the convergence of the SCSA for h — > 0. The 
work must be orientated at a second step to the comparison between the performance 
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of the SCSA and other signal analysis methods like Fourier transform or the wavelets 
and also to the generalization of the SCSA to other fields. 
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A Direct and inverse scattering transforms 

These appendices recall some known concepts on direct and inverse scattering 
transforms of a one-dimensional Schrodinger operator. For more details, the reader 
can refer to the large number of references on this subject for instance jT|, [3], [6], 

We consider here the spectral problem of a Schrodinger operator H\{— V), given 
by 

--X + V(x,M = k 2 ^, keC + , iGl, (46) 

dx z 

where the potential V such that V G B. For simplicity, we will omit the indice 1 of 

the spectral parameters in the following. 

For k 2 > 0, we introduce the solutions ip± of equation ( H6l) such that 



j, (kx) - I T(/C)e ^ X ~* ~°°' (47^ 

y-^'X) ~ \ e ~ ikx + R r (k)e +ikx x^+oo, K ' ' 

l (h \ - I T ( k ) e+lkx x -> +°°> (A ^ 

where T(k) is called the transmission coefficient and Ri^(k) are the reflection coeffi- 
cients from the left and the right respectively. The solution ^_ for example describes 
the scattering phenomenon for a wave e~ lkx of amplitude 1, sent from +oo. This 
wave hit an obstacle which is the potential so that a part of the wave is transmitted 
T{k)e~ lkx and the other part is reflected R r (k)e +lkx . ip + describes the scattering 
phenomenon for a wave e +thx sent from — oo. 

For k 2 < 0, the Schrodinger operator spectrum has N negative eigenvalues de- 
noted —k 2 , n =!,-•• ,N. The associated L 2 -normalized eigenfunctions are such 
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that 

ipn(x) = ci n e~ KnX , x ->■ +00, (49) 

iP n (x) = (-l) N ~ n c rn e +K - x , x ->-oo, (50) 

q„ and c rn are the normalizing constants from the left and the right respectively. 
The spectral analysis of the Schrodinger operator introduces two transforms: 

• The direct scattering transform (DST) which consists in determining the so 
called scattering data for a given potential. Let us denote Si(V) and S r (V) the 
scattering data from the left and the right respectively: 

Sj(V) := {Rj(k), K n , c jn , n = l,--- ,N}, j = l,r, (51) 

where j = r if j = I and j = I if j = r. 

• The inverse scattering transform (1ST) that aims at reconstructing a potential 
V using the scattering data. 

The scattering transforms have been proposed to solve some partial derivative 
equations for instance the KdV equation |10| . 

B Reflect ionless potentials 

Deift and Trubowitz [6] showed that when the Schrodinger operator potential V 
satisfies hypothesis ([2]) then it can be reconstructed using an explicit formula given 
by 

N 2i f + °° 

V(x) = -A s ^K n ^ 2 n (x) + — kR r(l) {k) fl(k,x)dk, 16I. (52) 

n =i n ^-°° 

This formula is called the Deift- Trubowitz trace formula. It is given by the sum of 

two terms: a sum of K n ip^ that characterizes the discrete spectrum, and an integral 

term that characterizes the contribution of the continuous spectrum. 

There is a special classe of potentials called reflectionless potentials for which the 

problem is simplified. A reflectionless potential is defined by Ri^(k) = 0, VA; 6 R. 

According to the Deift-Trubowitz formula, a reflectionless potential can be written 

using the discrete spectrum only, 

A' 



V(x,t) = -A^2K n ^{x,t), (53) 



n=l 



2(3 



C An infinite number of invariants 

There is an infinite number of time invariants for the KdV equation given by 
the conserved quantities |10| . |12| . |26j . Let us denote these invariants I m (V), m = 
0, 1, 2, • • • . They are of the form 

+1 2m + l f +ao n , xr dV d 2 V 



UV) = ("l) m+ '-^rj^ P m{V, ^, ^ ■ ■ ■ )dx, (54) 

where P m , m = 0, 1, 2, • • • are known polynomials in V and its successive derivatives 
with respect to x G M [3]. 

A general formula relates I m (V) to the scattering data of Hi(—V) [3], |12j . |26) 
as follows: 

N 9m 4- 1 r +OQ 

UV) = J2< m+1 + ^&- / (-k 2 ) m ^(l-\R ril) (k)\ 2 )dk, (55) 

n=l ^ J ~°° 

m = 0, 1, 2, • • • . So, for m = 0, P (V, • ■ ■ ) = V, we get with (153)1 and ([55]): 

/+oo -" 1 /"+oo 

V{x)dx = -A^K n -- I ln(l- |i? r(/) (fc)| 2 )^. (56) 

-oo n=1 7T y_oo 
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